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Abstract. We compute an the genus 1 correction to free energy of Hermitian two-matrix model 
in terms of theta-functions associated to spectral curve arising in large N limit. We discuss the 
relationship of this expression to isomonodromic tau-function, Bergmann tau-function on Hurwitz 
spaces, G-function of Frobenius manifolds and determinant of Laplacian in a singular metric over 
O \ spectral curve. 
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1 Two-matrix models: introduction 

O ■ 
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where the integral is taken over all independent entries of two hermitian matrices M\ and M2 such that 
^ . the eigenvalues of Mi are concentrated over a finite set of intervals (cuts) with given filling fractions. 
This integral is to be understood as a formal asymptotic series in N and in the coefficients of the 
two potentials V\ and V%. As a formal series, the questions of convergence of the matrix integral is 
irrelevant, and the model can be extended to matrices with eigenvalues constrained on contours in the 
complex plane. 

Such asymptotic series play an important role in physics, as generating functions of statistical 
physics on random discretized polygonal surfaces, i.e. a simplified model of euclidean 2D quantum 
gravity [31 E]. The large N expansion F = ^2q =0 N~ 2G F G (N is the matrix size), called topological 
expansion, is one of the cornerstones of the theory, since F G has the meaning of generating function for 
random discretized polygonal surfaces of genus G. Double scaling limits of these models correspond 
to statistical physics models on continuous surfaces, with conformal invariance properties. Matrix 
models thus provide realizations of minimal (p, q) conformal models. The 1-matrix model was shown 
to correspond to pure gravity (i.e. q = 2), and the 2-matrix model was introduced as it produces all 
(j), q) minimal models. 

Recently, the interest in large TV" matrix models was renewed as it was understood [27] , that the 
large N free energy of matrix models is the low energy effective action for some string theories. The 
computation of 1/N 2 expansion for both one-matrix and two-matrix models is based on the loop 
equations, which was first derived for 1-matrix 1-cut in [Sj, then for 1-matrix 2-cuts in [3E3) and 



In this paper we study the partition function of multi-cut two-matrix model: 



1 



recently derived in |12l 1131 II Ij for 1-matrix model multicut, and in (141 115j for two- matrix model 1 
and 2-cuts. Here, we will extend the results of [141 115j for an arbitrary number of cuts i.e. for an 
arbitrary (up to maximal) genus of the spectral curve. 
Writing down polynomials V\ and V 2 in the form 

di+l d 2 +l 
k=l k=l 

we shall use the following standard notations for operators of differentiation with respect to their 
coefficients: 

dl+l r d 2 +l 

- £ y^kd^ . (1.3) 



x kd u . 



SVi(x) ,■ ^ u * ' SV 2 (y) , , , 

These notations will be used below to shorten some of the formulas; by definition the equality 
8F 



5Vi{x) 



H(x) means that = — !— / x k H(x)dx , k = 1, . . . , d\ + 1 ; (1.4) 

du k 2mk J x=00 



a detailed discussion of this notation is contained in [^j. In fact, formally it is much more convenient 
not to cut the functions V\ and V 2 to polynomials, but instead consider the Laurent series 

00 00 

k=l k=l 



In this case we have the formal relations 

5Vi{x) 1 SVi'{x) 



8V\(x) x — x : 6V\(x) (x — x 
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(1.6) 



which are implicitly used in the derivation of loop equation. However, the convergency problem with 
considering all coefficients in the infinite sums ()1.5|) to be independent variables forces us to understand 
all relations involving the operators 5/5V\{x) and 6/5V2(y) in the sense of 1)1.4(1 . 
Consider the resolvents (also understood as formal power series) 

As a corollary of 1)1. 6|) , the free energy of two- matrix model ()1.1)) satisfies the following equations with 
respect to coefficients of polynomial V±: 

5F 5F 

W(s), TT7^T=W(y), (1-8) 



SVi(x) v ' SV 2 (y) 

valid in the sense of 1)1.4)1 . 

Assuming existence of 1/A^ 2 expansion, the equations 1)1.8)1 were solved in in the zeroth order 
in terms of holomorphic objects associated to the "spectral curve" which arises in iV — > 00 limit. The 
next coefficient F 1 was computed in )14j if the genus of spectral curve equals zero, and in ^3] if the 
genus equals one. The main result of this paper is an expression for F 1 for an arbitrary genus of 
"spectral curve", which we find using loop equations. We compute -F 1 using the algebro-geometric 
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framework (the spectral curve and corresponding machinery) which arises already in zeroth order 
approximation . 

The spectral curve is defined by the following equation: 

£°(x, y) := (V^x) - y)(V 2 '(y) - x) - V°(x, y) + l = 0, (1.9) 

where polynomial of two variables V°(x, y) is the zeroth order term in 1/N 2 expansion of polynomial 

1 ' U > N \ x - Mi y - M 2 I V ; 

the point P of this curve is the pair of complex numbers (x,y) satisfying (|1.9|) . 

The spectral curve (|1.9j) comes together with two meromorphic functions f{P) = x and fif(P) = y, 
which project it down to x and y-planes, respectively. These functions have poles only at two points 
of L, called ooj and oo g : at 00/ function f(P) has simple pole, and function g(P) - pole of order 
d\ with singular part equal to Vi(f(P)). At oo g the function g(P) has simple pole, and function 
f(P) - pole of order d 2 with singular part equal to V 2 '(g(P)). Therefore, one gets the moduli space 
A4 of triples (L,f,g), where functions / and g have this pole structure. The natural coordinates on 
this moduli space are coefficients of polynomials V\ and V 2 and g numbers, called "filling fractions" 
e a = §. gdf ', where a a are (chosen in some way) canonical cycles on L. 

Denote the zeros of differential df by Pi, ... , P mi (mi = d 2 + 2g + 1) (these points play the role of 
ramification points if we realize L as branched covering by function /(P)); their projections on /-plane 
are the branch points, which we denote we denote by h := f(Pj) ■ The zeros of the differential dg (the 
ramification points if we consider L as covering defined by function g(P)) we denote by Qi, . . . , Q m2 
(m 2 = d% + 2g + 1); their projections on g-plane (the branch points) we denote by \Xj := g(Qj)- We 
shall assume that our potentials V\ and V 2 are generic i.e. all zeros of differentials df and dg are simple 
and distinct. 

If is well-known ^fj] how to express all standard algebro- geometrical objects on L in terms of the 
previous data. In particular, the Bergmann bidifferential B(P,Q) = dpdQlnE(P,Q) (E(P,Q) is the 
prime-form), can be represented as follows: 

(see [16] for the proof). The Bergmann bidifferential has the following behaviour near diagonal P — ► Q: 

B{P ' Q) = { (t(P)-t(Q)) 2 + ^ Sb(P) + dT{P)dT{Q) ' (L12) 

where r(P) is some local coordinate; Sb(P) is the Bergmann projective connection (Sb(P) transforms 
as quadratic differential under Mobius transformations; under an arbitrary coordinate transformation 
an appropriate Schwarzian derivative is added to it). 

Consider also the four-differential D(P, Q) = dpdp lnP(P, Q), which has on the diagonal the pole 
of 4th degree: D(P,Q) = {6(r(P) - r(Q))~ 4 + 0(l)}dT{P){dT(Q)) 3 . From B(P,Q) and D(P,Q) it 
is easy to construct meromorphic normalized (all a-periods vanish) 1-forms on L with single pole; in 
particular, if the pole coincides with ramification point P&, the natural local parameter near P& is 
Xfc(P) = y/ f(P) — hi and the following objects: 

dx k (Q) Q=p k {dx k (Q)f Q=p k 
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are meromorphic normalized 1-forms on L with single pole at P k and the following singular parts: 
B(P, Pk) = {^py + ls B (P k ) + dx k (P) ; D(P, P k ) = + °(!)} dx k( p ) (1-14) 

as P — > Pfc, where SsiPk) is the Bergmann projective connection computed at the branch point P k 
with respect to the local parameter Xk{P)- 

Equations (jl.8|) in order l/N 2 look as follows (we write only equations with respect to V\): 

-Y\P) , (1.15) 



<^i(/(^)) 



where the Y 1 is the 1/iV 2 contribution to the resolvent W. The function Y 1 can be computed using 
the loop equations |14j . which leads to the following expression: 

The solution of l|1.15|) . (|1.16f) invariant with respect to the projection change (i.e. which satisfies also 
the required equations with respect to V2), looks as follows: 



mi 



F 1 = ± In \ rf{v d2+l ) 1 -^ dg(P k ) + **±* lnd 2 , (1.17) 



fc=l 



where tj is the so-called Bergmann tau-function on Hurwitz space, which satisfies the following system 
of equations with respect to the branch points l k '- 

JL lnTf = -L SB( P kh (1 .i8) 

The Bergmann tau-function (|1.18j) appears in many important problems: it coincides with isomon- 
odromic tau-function of Hurwitz Probenius manifolds ^Hl> an d gives the main contribution to G- 
function (solution of Getzler equation) of these Frobenius manifolds; it gives the most non-trivial term 
in isomonodromic tau-function of Riemann-Hilbert problem with quasi-permutation monodromies. 
Finally, its modulus square essentially coincides with determinants of Laplace operator in metrics 
with conical singularities over Riemann surfaces • The solution of the system (jl. 18f) was found in 
[2U] and can be described as follows. 

Define the divisor (df) = — 2ooj — (^2 + 1)009 + ^^1 Pk '■= Y^k=\ 2 r kDk- Choose some initial point 
P G C and introduce corresponding vector of Riemann constants K p and Abel map ^4a(Q) — Jp 
(w a form the basis of normalized holomorphic 1-forms on L). Since all zeros of differential df have 
multiplicity 1, we can always choose the fundamental cell £ of the universal covering of L in such a 
way that A((df)) = —2K P (for an arbitrary choice of fundamental domain these two vectors coincide 
only up to an integer combination of periods of holomorphic differentials), where the Abel map is 
computed along the path which does not intersect the boundary of C. 

The main ingredient of the Bergmann tau-function is the following holomorphic multivalued (1 — 
g)g/2-differential C(P) on L: 

t £S"'"' (p) -'"- (F) - (1 ' 19> 



4 



where 

W(P) := det^^gWw^iP)]] (1.20) 

denotes the Wronskian determinant of holomorphic differentials at point P. Introduce also the quantity 
Q defined by 

m+2 

Q = [df(P)]^C(P) n [E{P,D k )f^ , (1.21) 
fc=i 

which is independent of the point P £ L. Then the Bergmann tau-function (|1.18j) of Hurwitz space is 
given by the following expression: 

rre+n 

T/ = q2/3 tj [E(A,A)]^ ; (1-22) 

k,l=l k<l 

together with (jl,17|) this gives the answer for 1/N 2 correction in two- matrix model. 

If potential V<i is quadratic, integration with respect to M<i in (jl.ll) can be taken explicitly, and 
the free energy Q1.17JI gives rise to the free energy of one-matrix model. The spectral curve L in this 
case becomes hyperelliptic, and the formula (|1.17|) gives, using the expression for -ry obtained in |23| : 

1 f 23+2 1 

F 1 = - In { A* (det A) 12 ft g'(l k )\ , (1.23) 



k=l 



where Ik, k = 1, . . . , 2g + 2 are branch points of L; A is their Wronskian determinant; A is the matrix 
of a-periods of non-normalized holomorphic differentials on L. 

The paper is organised as follows. In section 2, following [T3|, we write down the loop equations 
for two-matrix model, and discuss the spectral curve and associated objects which arise in the zeroth 
order in 1/N 2 expansion. Here we derive also some new variational formulas, which will be used 
later in computation of 1/N 2 correction to free energy. In section 3 we solve the loop equations in 
1/N 2 approximation. Here we also express F 1 in terms of Bergmann tau-function on Hurwitz spaces 
introduced in |18[ I26|. In section 4 we recall the explicit expression for Bergmann tau-function [20], and 
find its transformation law under the change of projection of the spectral curve to CP 1 . This allows to 
get the formula for F 1 which satisfies the full set of variational equations with respect to polynomials 
V\ and V-2- In section 5 we derive variational equation of F 1 with respect to filling fractions. In 
section 6 we discuss the links between F 1 and other related objects: determinant of Laplace operator, 
G-function of Frobenius manifolds and isomonodromic tau-function of fuchsian system with quasi- 
permutation monodromies. 

2 Loop equations: leading term 

Introduce the function 

Y{x) = V({x) - W(x) (2.1) 
In terms of function Y equations (|1,8|) for free energy can be written as follows: 

SF 



as well as (|1.8|) . valid in the sense of ()1.4j) . 



Vi'(s) - Y(x) , (2.2) 
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To make use of variational formula (|2.2I) we need to get some information about the function Y(x). 
This information is in principle contained in the loop equations, which follow from reparametrization 
invariance of the partition function (see [E] f° r details). To write them down, apart from 

resolvent W{x) (|1.7j) . we need to introduce the following objects: 



Polynomial V(x,y): 



Vix y) - l/ ^l-^W^-W ) (2 3) 



• Polynomial £(x, y) 

£(x,y) :=(y 1 (x)-y)(y 2 (y)-x)-P(x,y) + l (2.4) 

• Function U(x,y), which is a polynomial in y: 

• Function £/(x, y, z), which is also a polynomial in y: 

Now we are in position to write down the loop equation 

U{x,y) = x-V 2 {y) + 77^-^72 ^FT\ 2 - 7 ) 

y — Y[x) N z y — Y(x) 

which arises as a corollary of reparametrization invariance of the matrix integral (jl.lj) |14j . 

The residue at y = 5^(ic) of (|2.7jl leads to the following loop equation (for polynomials of degree 3 
this equation was first derived in 5 ) for function Y{x) := Vi'(x) — W(x): 

£°(x,Y(x)) = j^U(x,Y(x),x) . (2.8) 

To use the loop equation effectively we need to consider the 1/N 2 expansion of all of their ingredients. 

2.1 Leading order term: algebro-geometric framework 

Assume that the function Y admits an expansion into a power series in 1/N 2 : 

Y(x) = Y° + j^Y 1 + ... . (2.9) 

Then in the leading order the master loop equation (|2.8|) turns into algebraic equation in two variables: 
x and Y(°\x): 

£(x,Y°(x)) = 0, 

where 

£°(x,y) = (li'(x) - y)(V 2 '(y) - x) - V\x,y) + 1 . (2.10) 
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The polynomial equation 

£°(x,y) = (2.11) 

defines an algebraic curve L of some genus g, which we call "spectral curve" (if the spectral curve 
is non-singular, it has "maximal genus" equal to d\d 2 — 1); the point P of this curve is a pair of 
complex numbers (x,y) satisfying the polynomial equation (|2.11|) . Therefore, Y° can be considered 
as multi- valued function of x. The curve L comes together with two meromorphic functions on it: 
function f(P) = x and function g(P) = y (= Y (x)). Since polynomial V (|2.<3|) and function £ (|2.4|) 
are symmetric with respect to substitution x «-> y, V\ «-> V2, the same algebraic curve appears if we 
write down the loop equations for X(y) : = V 2 '(y) — svf( y ) ■ 

Analytical properties of functions f(P) and g(P) on L are well-known (see [161 117j and references 
therein). Namely, f(P) and g(P) are meromorphic functions on L having poles only at two marked 
points 00 f and oo g with the following pole structure: function f(P) has simple pole at 00 j and pole 
of order d\ at oo g ; function g(P) has simple pole at oo g and pole of order d 2 at 00/. Therefore, near 
00 f we can write the singular part of g(P) as polynomial of f(P); near oo 5 we can represent the 
singular part of f(P) as polynomial of g{P); coefficients of these polynomials are given by V± and V 2 ', 
respectively: 

g(P)=V 1 '(f(P))-j± p - )+ 0(f- 2 (P)) as P^oo/, (2.12) 

f{P) = V 2 '{g{P))-^ F) +0(g- 2 {P)) as P -+ oc g . (2.13) 

The dimension of the moduli space of triples (L, /, g) satisfying these conditions equals d\+d 2 +g+2. 
Let us choose on L a canonical basis of cycles (a a , b a ). Then coordinates on the space M. can be chosen 
as follows: 

• d\ + 1 coefficients u\, . . . , u^+i of polynomial V\ . 

• d 2 + 1 coefficients v\, . . . , Vd 2 +i of polynomial V 2 ■ 



The "filling fractions" 



1 

2vri 



gdf . (2.14) 



do 



In strictly physical situation potentials V% and V 2 should be such that, considering L as a covering 
defined by function /, one can single out the "physical" sheet (which includes point 00 f) such that all 
a-cycles lie on this sheet and each a-cycle encircles exactly one branch cut (all corresponding branch 
points must be real if potentials V\ and V 2 are real). Similar requirement comes from g-projection 
of L. However, here we don't impose these "physical" requirements i.e. consider the "analytical 
continuation" of physical sector, in the spirit of |27j . 

Nevertheless, the sheet of the curve L (realized as d 2 + 1-sheeted branched covering by function 
/), which contains the point 00 j, is called the "physical" sheet; the physical sheet is well-defined at 
least in some neighbourhood of 00 f. Fixing some splitting of L into d 2 + 1 sheets, we denote by x^ 
(k = 1, . . . , d 2 + 1) the point of L belonging to kth sheet such that f(x^) = x; we assume that point 
x^ belongs to the physical sheet of L i.e. x^ — > 00/ as x — > 00. 

The polynomial £°(x,y) defining the spectral curve L (|2.11j) can also be rewritten as follows: 

d 2 +l 

£°(x,y) = -v d2+1 J] {v-g{x {k) )) (2-15) 

k=l 
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The proof of (|2.15|) is simple: function £ is given by l|2.1U|) : since V° is a polynomial of degree d% — 1 
with respect to y, function £ is a polynomial of degree d 2 + 1 in y; its zeros are Y°(x( fc )) by definition 
of points x( k \ Comparison of coefficient in front of y d2+1 leads to (|2.15|) . 

2.2 Some variational formulas 

If a Riemann surface is realized as a branched covering of Riemann sphere, the branch points can be 
used as natural parameters on the moduli space, and it is easy to differentiate all objects introduced 
above with respect to the branch points. The answer is given by Rauch variational formulas (j7j, 
for a simple proof see [23 )• However, on our moduli space the set of natural coordinates is given 
by coefficients of polynomials V\ and V 2 and filling fractions. To differentiate all interesting objects 
with respect to these coordinated we need to know the matrix of derivatives of branch points (we 
shall consider only {Ik}) with respect to coefficients of V{, V 2 and filling fractions. This matrix was 
computed in [EJ; below we rederive some of these formulas, and prove new variational formulas, 
required in our context. 

In ^S] equations (|2.2|) . together with their counterpart with respect to V 2 (y), were solved in the 
leading term i.e. it was found the solution of the system 

5F° 

= Vt'ifiP)) - g(P) 



5F° 



= V 2 '(g(P))-f(P) 

9{P) 



6V 2 (g(P))\ 

which a posteriori turns out to satisfy also the following equations with respect to filling fractions: 

OF 



^_ = T a :=j> g{P)df{P) . 



To find solution of the equations (|2.2jl in order 1/iV 2 , together with their counterpart with respect 
to V2(y)) we shall need 

Lemma 1 The following variational formulas take place: 



sviUiP)) 



g'(P k )df(P)=B(P,P k ) , (2.16) 



S{f/{Pk)} df(P) = \ (D(P,P k ) - 9 -^B(pP k ) } (2.17) 



8V 1 (f(P)) f(P) Jy ' 4 1 v ' K > g'(P k ) 
Proof. We start from formula (|1.11|) for the Bergmann bidifferential: 



B{p - Q)= ^mkL> dmd,(Q) - (2 ' 18> 



We want to rewrite this formula in the limit Q — > P k using the local parameter x k (Q) = \J f(Q) — h- 
As the first step we notice that for any coordinate t on our moduli space we have the following identity: 

9t(Q)\f(Q)df(Q) = gt(Q)\ Xk ( Q )df(Q) - ft(Q)\ Xk{ Q) d 9(Q) > ( 2 - 19 ) 
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which follows from differentiation of composite function g(t, f(x k ,t)) with respect to t using the chain 
rule. In particular, 

5g{Q) (Q)df(Q) = J^L ,MQ)~J^L ,MQ)- (2-20) 



5V x {f{P)) f^' J SVi(f(P)) SViiHP)) 



x k (Q) 



Consider now first several terms of local expansion of g(Q), dg{Q) and B(P,Q) as Q — » P k (prime 
denotes derivative with respect to x k := x k (Q)): 

g(Q) = g(Pk)+g'(P k )x k + ... , (2.21) 

MQ) = W(Pk) + g"(P k )x k + \g'"{P k )xl + ... }dx k , (2.22) 

B(P, Q) = {B(P, P k ) + B'{P, P k )x k + \b"{P, P k )x\ + ... }dx k . (2.23) 

Taking into account that f(Q) = x\ + l k , and substituting these relations into (|2.2U|) . we get in the 
order zero the formula (|2.16|) . 

The first order terms give relation which defines the dependence of g{P k ) on {u k }: 

2 w§k ~ mrn/ w ) di(p) - B ' {p - Pt) ' {2M) 



we present this relation only for completeness, since it will not be used below. 



Finally, collecting the coefficients in front of x|, we get 



2 tg'(Pk) _ 1 si k gl „ _ 1 B"(p,p k ) 



SViifiP)) 2 5V 1 (f(P)f y K > 2 df(P) ' 
which leads to IJ2.17J) after using (|2.1tijl . 

3 Solution of loop equation in 1/N 2 approximation 

The main goal of this paper is to find function F 1 on our moduli space which satisfies the equation 

8F 1 



SVxix) 



-Y\x) , (3.1) 



where Y 1 (x) should be determined from l/N 2 expansion of the loop equation (|2.8|) . The equation 
()3.1|) is valid in a neighbourhood of the point oo f i.e in a neighbourhood of the point x = oo on the 
"physical" (with respect to variable x) sheet of the spectral curve L. The same function should 
satisfy the equation 

ftF 1 

S45) -*'<»>■ ^ 

where function X l (y) should be found from writing down the loop equation with respect to matrix 
M2 in a neighbourhood of point oo 9 . We shall first solve ()3.1|) . and then check the symmetry of the 
expression with respect to the change of projection f g. 
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To express Y 1 in terms of the objects associated to the spectral curve L we consider the 1/iV 2 
term of the master loop equation (|2.8[) : we have: 

£[x, Y(x)) = £°(f(P),g(P) + ^Y\P) + ...) + -^S\f(P),g(P)) + . . . , (3.3) 

as P — > oof, where, as before, in a neighbourhood of ooj, f(P) = x; g(P) = Y°(x). The 1/N 2 
expansion of £(x, y) looks as follows: 

S(x, y) = £°(x, y) + ±£ \x, y) + . . . ; (3.4) 

since £ 1 (x,y) = —V 1 (x,y), we can further rewrite this expression in a neighbourhood of point oof as 
follows: 

S(x,Y(x)) = S°(f(P),g(P) + ±{£\f(P),g(P)) + Y\P)£$(f(P), g(P))} + . . . . (3.5) 

Therefore, the 1/iV 2 term of master loop equation (|2.8j) gives 

U°(f(P),g(P)J(P)) = £\f(P),g(P)) +Y 1 (P)£° y (f(P),g(P)) 

as P — > oo/ , or 

To make this formula more explicit we need to express U°(f(P),g(P), f(P)) in terms of known objects 
using the loop equation ()2.7j) . According to definition of U°(x,y, z) we have: 

«o (w ) = -^# (3 . 7) 
dVi{z) 

On the other hand, the zeroth order term of ()2.7|) gives: 

U°(x,y)=x-V 2 '(y) + £ ° {X ( ,y L (3.8) 

y-g[x^>) 



(as before, x^ denotes a point on the physical sheet of L). Therefore, 

(3.9) 

Using the form (|2.15j) of the polynomial £°(x,y), we can further rewrite this expression as follows: 



j*, x 6£°(x,y)/5V 1 (z) £°(x,y) Sgjx^) 
[X,V,Z) ~ v-g(xW) (y-g(xU)r SVi{z) 



*W)_ 1 . (3 ,0) 



W*) y-5(^ (fc) ) ' 

Substituting this expression into (|2.15jl . we get 



i ( , u , 



fc=2 
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Substituting z = x = f(P) and taking the limit y — > (/(a^ 1 )), we have: 

U°(f(P),g(P),f(P))=£° y (f(P),g(P)) ]T g X — _ (3.12) 

as P = — > ooj. Now Q3.6)) can be rewritten as follows: 

yl f m = g(/(W) V W) 1 

lJ £§(f(P),g(P)) + Q ^ p ^ )=m 5V 1 (f(P))g(P)-g(Q) { ' ' 

as P — > oof, which can be further transformed, using the formula (|l.llj) for the Bergmann bidifferen- 
tial: 

„,,„,„ WW5)„ PU V- B(P,Q) 1 

Y imiF)= wmm^ lF)+ Q ^^ tm w) 9(P)- S (Q) : (3 ' 14) 

now we see that the 1-form Y 1 (P)df(P) can be analytically continued from a neighbourhood of ooj 
to the whole L. 



Lemma 2 Let the spectral curve L \2.11\) be non-singular. Then the 1-form Y 1 (P)df(P) {3. 14\ ) is a 
meromorphic 1-form on the spectral curve L with poles up to fourth order only at the branch points 
Pk i.e. at the zeros of differential df{P). 



Proof. Let us verify the non-singularity of the first term 

V 1 {f{P),g{P)) 



S°(f(P),g(P)) 



df(P) , (3.15) 



of (|3.14|) everywhere on L. For finite f(P) the 1-form Q3.15[) can be singular only at the zeros of 
£y(f(P),g(P)), which, if the curve L is non-singular, are by definition the branch points Pk] these 
zeros are of the 1st order and are canceled by the zeros of df(P) at the branch points. 

To study behaviour of (|3.15j) at ooj and oo g we mention that the polynomial V(x,y) (|2.3jl (and, 
therefore, also its first correction P 1 (x, y)) is of degree d± — l with respect to x and d 2 — 1 with respect 
to y. However, we can say a bit more about V 1 (x,y). Namely, the coefficient of V(x,y) in front of 
yd.2—1 e q Ua ls u c i 1 +iVd 2 +i, which does not have any higher corrections. Therefore, the coefficient 



x 



di—l„.d,2 



of V 1 (x,y) in front of x dl l y d2 1 vanishes. 

Now consider the behaviour of the 1-form ()3.15|) near oof. We have 

£° y (f(P),g(P)) = -(V 2 '{g(P)) - f(P)) - W{f(P)) - g{P))V 2 "{g(P)) - V ° y (f(P),g(P)) , 



which has pole of order did 2 near oof as corollary of asymptotics p. 12)1 of function g(P) near oof. The 
1-form df(P) has at oof the pole of second order. The most singular contribution by V 1 (f(P), g(P)) 
at oof comes from the term f dl ~ 2 (P)g d2 ~ l (P); it has the pole of order d\ — 2 + d\(d 2 — 1) = d\d 2 — 2. 
Summing up, we see that (j3.15j) is non-singular near oof. 

Consider the 1-form (|3.15|) near oo g . At this point df{P) has pole of order d 2 + 1; the main contri- 
bution to £°(/(P),5(P)) is given by the term (Vi{f{P)) - g(P))V 2 " (g(P)), which has pole of order 
^1^2+^2 — 1. Finally, the main contribution to V 1 (f (P) , g{P) ) comes from the term g d2 " 2 (P)/ dl - 1 (P), 
which has pole of order d\d 2 — 2. In total 1)3. 15j) is non-singular at oo g , too. 
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Consider now the second term of (|3.14j) : 

B(P,Q) 1 



df(Q) 'g(P)-g(Q) 
Q¥=P ■ f(Q)=f(P) 



(3.16) 



The Bergmann bidifferential is singular (has second order poles) only at coinciding arguments, i.e. 
when P coincides with one of the branch points Pk- The denominator g(P) — g(Q) vanishes only if P 
coincides with Q, (i.e. both of them coincide with one of the branch points Pk). It is slightly more 
complicated to see that zeros of df(Q) don't give any poles outside of Pk- Obviously, df(Q) is singular 
if P — > P k and Q = P*, where P* is another point such that f(P*) = f(P) and P* — > P k as P — ► Pk- 
However, df(Q) is also singular if Q coincides with one of the branch points Pk, while P remains on 
some other sheet, and does not tend to Pk as Q —* Pk- In this case in the sum ()3.16j) we have two 
singular terms (with poles of first order), which correspond to Q and Q*; however, the residues of 
these terms just differ by sign, and, therefore, the total sum (|3.16j) remains finite outside the branch 
points Pk and infinities oof and oo g . 

As P — » ooj, all points Q in (|3.16|) tend to oo g ; thus all df(Q) have pole of order c?2 + 2; all other 
terms remain non-singular and non- vanishing. Therefore, (|3. 16|) has zero of order c?2 + 1 at oo/- 

As P — > oo 9 , the situation is slightly more complicated. If we enumerate the sheets of L such that, 
as x — > oo, x^ — > oof, and x^ 2 \ . . . ,x^ d2+1 ^ — ► oo 5 , and assume that P = x^ 2+1 \ then H3.16j) can be 
split as follows: 

B( x (i) ja; («fa+i)) i * fi(i(i),x(*+i)) 1 

As x — > oo, the first term in ()3.17j) has zero of order two (df(x^) has pole of order two, other 
multipliers remain non-singular and non- vanishing) . In each term of the sum in (|3.17j) the Bergmann 
bidifferential has pole of second order as x — » oo; however, df{x^) has pole of order c?2 + 1 ; and 
g{x <yd2+l ' > ) — g(xV') has simple pole as x — > oo; therefore, the whole expression (|3.17|) is non-singular 
(and even vanishing) as x — > oo. 

Remark 1 The condition of non-singularity of the spectral curve (j2.11jl made in lemma means 
in physical language that the spectral curve has maximal possible genus equal to d\d2 — 1 for given 
degrees of polynomials V± and If the genus of the spectral curve is less than the maximal genus, 
it must be singular; then the non-singularity of 1-form Y 1 (P)df (P) at the double points can not be 
verified rigorously. However, this non-singularity is suggested by physical consideration: since we 
assume that at the double points one does not have any eigenvalues of M\ or M2 in large limit 
(i.e. corresponding filling fractions are equal to zero), there is no physical reason for corresponding 
resolvents to be singular at these points in large N limit. Therefore, in the sequel we shall assume 
that Y 1 (P)df (P) is non-singular outside of branch points of L both for maximal and non-maximal 
genus. We should mention that this assumption was also made (explicitly or implicitly) in the previous 
papers EJ E3 Eg . 

The singular parts of Y l (P)df(P) at the branch points Pk can be found from (|3.14j) . If, say, 
P — > Pk, the only term in 1)3. 14J) which contributes to singular part at Pk corresponds to Q = P*. 
Thus 

r ' ww= w?ffiV) +Wi " P ~ Pk - (3 ' 18> 
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Consider the local expansion of all ingredients of this expression as P — > P k in terms of the local 
parameter x k (P) = \J 'f(P) - l k : 

g{P) = g{P k ) + x k {P)g'{P k ) + \x 2 k {P)g"{P k ) + \xl{P)g'" {P k ) + ... , 

2 b 

g(P*) = g(P k ) - x k {P)g'(P k ) + l -x 2 k (P)g" {P k ) - \x\W \P k ) + . . . , 

df(P*) = 2x k (P)dx k (P) , 

B{P ' n = ( (2x k (P)f + l SB{Pk) + • • • ) dx k (P)(-dx k (P)) . 

i i ( 1 _x^ g _^p k y 



Therefore, 



g{P)-g{P*) 2x k (P)g'(P k ) \ 6 g'(P k 

and, as P — > P k , 

b(p,p*) i \ i s B \ i nMl ,, , p , 



df(P*) g(P)-g(P*) \ 16xl(P)g'(P k ) \96g»(P k ) 24g'(P k )J x 2 k (P) 

(3.19) 

Since, according to our assumption, the 1-form Y 1 (P)df(P) is non-singular on L outside of the 
branch points, we can express this 1-form in terms of differentials B(P,P k ) and D(P,P k ) (|1.13|) using 
their behaviour near P k \ 

r(l, ™ ) -|{-^ D(p - Pt)+ [™-wK)] B(p - p 4 ; (32o) 



as a result we rewrite the equation (j3.1j) for F 1 as follows: 

"" r ' g'"{P k ) ; Snip,, 



5F 1 mi ( i 



Pk) + 



+ 



96 5 ' 2 (P fc ) 2Ag'(P k ) 



B(P,P k )\ . (3.21) 



Proposition 1 A general solutions F 1 of the system can be written as follows: 

1 I ( mi ) 

F 1 = ~]nT f + -]nUJg , (P h )\+C({v k },{e a }) (3.22) 

where C({v k }, {e Q }) is a function on our moduli space depending only on coefficients of polynomial V2 
and filling fractions {e Q }; function Tf (the Bergmann tau-function on Hurwitz space) is defined by the 
system of equations with respect to the branch points {l k }: 

|-hr / = -ifl,(J Wi (3.23) 

function Tf depends on coordinates {u k , v k , e a } as a composite function. 

Proof. The derivative of tj with respect to Vj.(/(P)) is computed by chain rule using variational 
formula (|2.16|) : derivatives of g'(P k ) with respect to V\{f{P)) are given by (|2.17|) . Collecting all these 
terms together we see that derivative of (j3.22|) coincides with (|3.21|) . 

Therefore, to compute -F 1 it remains to find the Bergmann tau-function Tf and to make sure that 
"constant" C({v k },{e a }) is chosen such that the final expression is symmetric with respect to the 
change of "projection" i.e that F 1 satisfies also equations (|3.2j) . 
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4 and Bergmann tau-function on Hurwitz spaces 

4.1 Bergmann tau-function on Hurwitz spaces 

Here, following ^2U], we discuss the Bergmann tau-function on Hurwitz spaces for the stratum of 
Hurwitz space arising in the application to two-matrix model. 

The Hurwitz space H g ^ is the space of equivalence classes of pairs (L, /), where L is a compact 
Riemann surface of genus g and / is a meromorphic functions of degree N. The Hurwitz space is 
stratified according to multiplicities of poles of function /. By H g ^(k\, . . . , k n ), where k\ + ■ • • + k n = 
N, we denote the stratum of H g ^ consisting of meromorphic functions which have n poles on L with 
multiplicities k\, . . . , k n . (In applications to two-matrix model we need to study the tau-function on 
the stratum H g ^(l, N — 1), on which the function / has only two poles: one simple pole and one pole 
of order N — 1.) 

Suppose that all critical points of the function / are simple; denote them by Pi,...Pm (to = 
2N + 2g — 2 according to Riemann-Hurwitz formula); the critical values l k = n(Pk) can be used as 
(local) coordinates on H g at(&i, . . . , k n ). Function / defines the realization of the Riemann surface 
L as an iV-sheeted branched covering of CP 1 with ramification points Pi, ... , P m and branch points 
h = f(Pk)', the points at infinity we denote by ooi, . . . , oo n . In a neighbourhood of the ramification 
point Pfc the local coordinate is chosen to be x k '■= \/l — lk: k = 1, . . . , to; in a neighbourhood of the 
point ooj the local parameter is x m+ j := l~ l l k i 

The Bergmann bidifferential B(P,Q) has the second order pole as Q — > P with the local be- 
haviour Ot : B(P,Q)/{dx(P)dx(Q)} = (x(P) - x(Q))~ 2 + ±S B (x(P)) + o(l), where x(P) is a local 
coordinate; Sb(x(P)) is the Bergmann projective connection. 

We define the Bergmann r-function Tf(l\, . . . , l m ) locally by the system of equations Ij3.23j) : 

d 1 

— In Tf = -—S B {xk)\x k =o , k = l,...,m. (4.1) 

compatibility of this system is a simple corollary of Rauch variational formulas |21| . 

Consider the divisor of the differential df: (df) = Y^k=i r kDk where Dk ■= Pk , »"fc := 1 
for k = 1, . . . ,m and D m+ j = ooj , r m+ j = —(kj + 1) for j = 1, . . . ,n; the corresponding local 
parameters Xk, k = 1, . . . ,m + n were introduced above. 

Here and below, if an argument of a differential coincides with a point Dj of divisor (df), we 
evaluate this differential at this point with respect to local parameter Xj. In particular, for the prime 
form we shall use the following conventions: 

E(D k ,D l ) := E(P, Q)Vdx k (P)^d Xl (Q)\ P=Dkt Q=Dl , (4.2) 

for k,l = 1, . . . ,m + N . The next notation correspond to prime-forms, evaluated at points of divisor 
(df) with respect to only one argument: 

E(P,D l ) := E(P,Q)yJd Xl (Q)\ Q=Dl , (4.3) 

/ = 1, . . . , to + n; in contrast to E(Dk, Di), which are just scalars, E(P,Di) are — 1/2-forms with 
respect to P. 

Denote by w\, . . . ,w g normalized (<f wp = 5 a p ) holomorphic differentials on L; = <f> b wp is 
the corresponding matrix of 6-periods; 0(z|B) is the theta-function. 

Choose some initial point P £ £ and introduce corresponding vector of Riemann constants K p 
and Abel map *A a (Q) — Jp w a . Since zeros of differential df have multiplicity 1, we can always choose 
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the fundamental cell C of the universal covering of L in such a way that A{{df)) = —2K P (for an 
arbitrary choice of fundamental domain these two vectors coincide only up to an integer combination 
of periods of holomorphic differentials), where the Abel map is computed along the path which does 
not intersect the boundary of L. 

The key entry of the Bergmann tau- function is the following holomorphic multivalued (1 — g)g/2- 
differential C(P) on L: 

' ai,...,a g =l 3 

where 

W(P) := det^^gWw^ (P)\\ (4.5) 

denotes the Wronskian determinant of holomorphic differentials at point P. 
The following theorem is a slight modification of the theorem proved in [2U] . 

Theorem 1 The Bergmann tau-function j4-l\) °f Hurwitz space H g ^(k\, . . . , k n ) is given by the fol- 
lowing expression: 

m+n 

Tf = Q 2/3 I] [EiP k MfP (4.6) 

k,l=l k<l 

where the quantity Q defined by 

m+N 

Q = [df(P)]^C(P) [J [E{P,D k )]^¥± ; (4.7) 

k=l 

is independent of the point P £ L. 

The proof of the theorem is very similar to [2U|. The only technical difference is the appearance 
of higher order poles of function /. 

4.2 Dependence of Bergmann tau-function on the choice of the projection 

Theorem 2 Let Tf and r g be Bergmann tau-functions \4-b\ ) corresponding to divisors (df) and (dg), 
respectively. Then 

aj^^a® (4 . 8) 

where 

C = (4.9) 
is a constant independent of moduli parameters. 

Proof. As above, we assume that the fundamental cell C is chosen in such a way that A((df)) = 
A{{dg)) = —2K P . Denote divisors (df) and (dg) as follows: 

mi mi+2 

(df) = Y^P k - 2oo/ - (cfc + l)oo fl := £ r k D k , (4.10) 

fc=l k=l 
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rri2 m.2+2 

(dg) = Y^Q k -2oo g - (di + 1)00/ := ^ s k G k . (4.11) 

k=l k=l 

Since deg(df) = deg(tfg) = 2g — 2, we have Y^k=\ 2 r k = Y^k=t 2 s k = 2g — 2. Then, according to the 
expression ([4.7ft for the Bergmann tau-function, we have 

mi+2 mi+2 

(r f ) 12 = C 8 (P)[df(P)] 4 ^ 4 J] {E(D k ,D j )} 2r ^ H {^(^^)} rfc(4 " 49) , (4-12) 

k,j=l k=l 

where the values of prime- forms at the points of divisor (df) are evaluated in the system of local 
parameters defined by function / i.e. near P k the local parameter is x k = w f(P) — l k ; near 00/ the 
local parameter is x mi+ \ = 1//(P), and near 00 9 it is x mi +2 = [f (P)]~ l I d2 . 
Similarly, we have 

m-2+2 m 2 +2 

(r g ) 12 = C 8 (P)[dg(P)]^' 4 [] {E(G k ,G j )} 2s ^ f[ {E(P,G k )} s ^-^ , (4.13) 

k,j=l k=l 

where the values of prime-forms at the points of divisor (dg) should be evaluated in the system of 
local parameters defined by function g i.e. near Q k the local parameter is y k = y/g(P) — \x k ] near 00 f 
the local parameter is y m2 +i = l/d(P)} an d near oo g it is y m2 +2 = [diP)]" 1 ^ 2 ■ 
Therefore, 

Tf\ 12 = nK?w.^)} Vi f df(p) wiit 2 {E(p,G k )y^ 9 - 4 

T J W^H{E{Gk,G 3 )} 2 ^ \dg{P) Uklf{E(P,D k )y k J 

Using independence of this expression of the choice of point P we can split the (Ag — 4)th power in 
this formula into the product over points of divisor (df) + (dg) (whose degree equals exactly 4g — 4!). 
The subtlety which arises is that, evaluating the prime-forms and differentials df and dg at the points 
D k and G k we fix the local parameters (these local parameters at the points of (df) are defined via 
function /, and at the points of (dg) via function g as explained above). Since divisors (df) and (dg) 
have common points (00 f and 00 g ), in a neighbourhood of each of these points we introduce two 
essentially different local parameters, and it is important to remember in each case in which local 
parameter the prime-forms are computed. 

Another subtlety is that, being considered as functions of P, different multipliers in Q4.14JI either 
vanish or become singular if P € (df) + (dg); cancellation of these singularities should be accurately 
traced down. 

Consider the first "half" of this expression, namely, the product over P £ (df): 
f df(P) UT=t 2 {E(P, G k )} Sk \ 2g ~ 2 m f^ 2 \ df(P) YlT=f{E(^Gk)} 



TT run ^^ii^.^j (41B) 



\dg(p) UT=t 2 {E(P,D k )y k \ " l\ p^ldgiP) UT=t 2 {E(P,D k )y 

if {E(D h D k ) } -^ IT 2 ! hm f «\ V k nl Hm 
}AJ L\ {P^D k {E(P,D k )y k j 11 )P-+ Dl dg(P) 



(4.16) 
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The first product looks nice since it cancels out against the first product in the numerator of (|4.14l) . 
Let us evaluate other ingredients of this expression. We have = P k , = 1 for k = 1, . . . , mi, 
Arn+i = 00/ , k mi+1 = -2, D mi+2 = oo g , k mi+2 = ~{d 2 + 1). Therefore, 

Tr 2 / v V k 



lim {df(P)E 2 (P,D mi+1 )}\ l \ lim {d/(P)^ +1 (P, J mi+2 )} 1 ' ' f[ lim f (P) , 

(4-17) 

where we don't write oof and oo 9 instead of D mi+ i and -D mi +2, respectively, to remember that we 
need to use the system of local parameters related to f(P). The last term in (|4.17f) product is the 
easiest one: 

lim P d ! {P) „ = hm ^ = 2 . (4.18) 
p^p fe {P(P,P fc )} x k (P)^o x k 

In a similar way we evaluate the first term: 



lim {df(P)E 2 (P,D mi+1 )} = -l, (4.19) 

f- 

and the second one: 



P^P mi+ i 



n lim {df(P)E d * +1 (P,D mi+2 )} = -d 2 . (4.20) 

P^P mi +2 

It remains to evaluate the third product in (|4,16|) : 



nsk 



Z=l I yV ; J \2=1 / \all such that D t ^G k 

xf lim {P(P,G m2+2 )} dl+1 d 5 (P)) ( lim {P(P,G m2+1 )} 2 ^(P)V + . (4.21) 

Consider the first limit in (|4.21l) : 
Lemma 3 

lim f{P(P,G m2+2 )} dl+1 ^(P)) 2 = (dDiu^) 1 -^ (4.22) 



Proof. Two different local parameters at the point oo/ = D mi+ i = G m2 +2 which we need to use are 
x mi+1 (P) = f-\P) and y m2+2 (P) = g^'^iP). We have 



cv pr , x {Vm 2 +2{P) + ■ ■ ■) dx mi+1 (y m2+2 (P) + . . .) 

E(P, G m2+2 ) = — —===— = t (oo/j — . (4.23) 

d^/y m2+2 {P) V ay m2 +2 y'dx mi+ i(P) 



We have also g(P) = Vm 2 \ 2 \ thus 



d 5 (P) = -(diXwsr^ 1 (^^(00/)) dx mi+1 (P) . (4.24) 
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Taking in Q4.22|) the limit P — > D mi+ i we indicate that all differentials in the bracket should be eval- 
uated with respect to the local parameter x mi +i- Therefore, in (|4,22|) we ignore all factors dx mi+ \(P) 
and (|4,22|) turns out to be equal to 

(4) (^^(pof)) l ^ = (dl){u dl+1 f-^ , (4.25) 
\dx mi+ i J 

where we take into account that, as P —* ccj, g = u^+ix^ 1 + ...; thus (dy m2+ 2/dx mi+ i)(oo f) = 
(u dl+1 )- 1 ^. 

Consider now the second limit in l|4.21|) : 
Lemma 4 

lim {E(P,G m2+1 )} 2 dg(P) = -l (4.26) 

Proof. In analogy to ()4.22|) we have to evaluate the prime-form and the differential dg with respect to 
the local parameter related to function / i.e. with respect to x mi+ 2 

(P) = (/(P))- 1 ^ 2 , while the local 
parameter arising from function g is y m2+ i(P) = (g(.P)) _1 . We have near D mi+ 2'- 

h(P, G m2+1 ) - . — {oog) — . (4.27) 

VVma+H-P) V dVmz+l yj dx mi+2 (P) 

and 

dg(P) = d ( —±j=r) = -^±l(oo 9 ) %^ . (4.28) 

As before, substituting these expressions to (|4,2fi|) and ignoring the arising power of dx mi +2(P), we 
see that this limit equals —1. 

Substituting this —1, together with the answer (|4.25|) for the limit (|4.22|) . into (|4.21|) . and collecting 
all terms in ()4.16j) . we get 

\dg(P) Y\Ilt 2 {E(p,D k )y 

ris k 



Now, computing the second "half of (|4.14() i.e. taking the product analogous to (|4. 16f) over points 
of divisor (dg) , and taking the product with ()4.29|) , we get the statement of theorem |2J 

4.3 Bergmann tau-function and 

Theorem 3 The F 1 solution of equations EO)) . JUJ) is given by any of the following two 

equivalent formulas: 

if i mi 1^-1-3 

= ^hJrjWi) 1 ^ ]Jdg(P k )\ + ^±-\nd 2 + C (4.30) 
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= 24 ln p^ 1 ) 1 " 5 Udf(Qk)\ + -^p lndi + C . (4.31) 

Here Tf and r g are Bergmann tau-function \4-b\) built from divisors (df) and (dg), respectively; C is 
a constant. 

Proof. Prom the formulas (|4.8|) . (|4.9|) it follows that expressions (|4,3U|) and (|4.31|) define the same 
function. According to Proposition^ expression (|4.3U|) satisfies equations (|3.1j) . (j3.2Uf) with respect 
to coefficients of V\. Similarly, expression (j4.31j) satisfies analogous system (j3.2j) with respect to 
coefhcients of Vi . 

Remark 2 (higher order branch points) If potentials V\ and V2 are non-generic i.e. some (or all) of 
the branch points have multiplicity higher than 1, formula (|4.31j) should be only slightly modified. 
Namely, the expression for Bergmann tau-function (|4.6|) formally remains the same in terms of divisor 
of differential df (the zeros of df can now have arbitrary multiplicities). The expression for F 1 then 
looks as follows: 

F 1 = ^Injrf (W^nHp^} + ^lnrf 2 + C. (4.32) 

The proof of (|4.32|) is slightly more involved technically than the generic case and will be published 
separately. 



5 Equations with respect to filling fractions 

it is well-known (see for example \17\ ) that the normalized (§ a wp = 5 a b) holomorphic differentials 
can be expressed as follows in terms of / and g: 



2m Wa (P) ~ d9{P) 



(Sketch of the proof: differentiating (|2.14j) with respect to e/3, we verify the normalization conditions for 
differentials (|5.1|) . The 1-form gdf is singular at oof and oo g ; at 00/ we have g = Vi(f) — 1// + this 
singularity disappear since coefficients of V\ and V2 are independent of filling fractions. Singularities 
at branch points Pk of derivative of g with respect to e a get cancelled by zeros of df at these points. 
At 00 g we have: x = V2(g) — 1/ g + due to thermodynamic identity 



dg 



At 9f 

f de a 



dg . 



Since coefficients of V% are independent of e a , singularity of gdf at 00 g also disappears after differen- 
tiation.) 

To obtain equations for derivatives of with respect to the filling fractions we shall prove the 
following analog of lemma Q 

Lemma 5 The following deformation equations with respect to filling fractions take place: 

o e Jk = -2m— n — T (5.2) 
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d{g\P k )} Tri f g"'(P fc ) rpU 



de a 2 1 av KJ g'{P k ) 

Proof is parallel to the proof of l|2.16|) and (|2.17|) : from ()5.1() we have 

(><I(P) d f(P)-^l dg(P) = 2niw a (P). (5.4) 



de a 



xk(P) de a 



Substituting in (|5.4j) the local expansions (|2.21|) of g(P) and ()2.22|) of dg(P), and expansion of w a (P) 

w a (P) = K(P fe ) + w' a (P k )x k + ^p^-x 2 k + ... )dx k , (5.5) 
we get, since f(P) = x 2 k {P) + l k and df(P) = 2x k {P)dx k {P): 

(d ea g(P k )+x k d e J(P k ) + \d e J'{P k ) + ... )2x k dx k - d e J k (g'(P k ) + g"{P k )x k + ^g'"(P k )x 2 k + ... )dx k 

= 2m(w a (P k ) + w' a {P k )x k + ^Wa(P k )xl)dx k . 

The zeroth order term gives (j5.2|) . Collecting the coefficients in front of x\, and using H5.2() . we get 
(tO) . 

Theorem 4 Derivatives of function F 1 \4-3U\) , i4-31\) with respect to the filling fractions look as 
follows: 

^- = -1 Y\P)df(P) , (5.6) 

ae a Jb a 

where Y l df is defined by \3.2(J\) . 

Proof. The vectors of 6-periods of these 1-forms B(P, P k ) and D(P, P k ) can be expressed in terms of 
holomorphic differentials via the following standard formulas: 

I B(P, P k ) = 2niw a (Pk) , <f D(P, P k ) = 2™<(P fc ) . (5.7) 

Jb a Jb a 

Therefore, the 6-periods of the 1-form -Y (l \P)df(P) (l3~2"Ul) are given by the following expression: 

mi ' w'i{P k ) , g"'(P k )w a (P k ) | S B (P k )w a (P k ) } ^ 



- 1 Y^(p )d f(p) = ™£i-<™ + 9 "' {p i w ; [ 



k) 24g'(P k 



On the other hand, derivatives of F 1 ()4.3U|) with respect to e a can be computed using ()5.2|) . (|5.3|) and 
equations for Bergmann tau-function (j3.23|) . which also leads to (|5.8|) . 

6 F 1 of two-matrix model, isomonodromic tau-function, G-function 
of Frobenius manifolds, and determinant of Laplace operator 

Here we outline some links between the expression (j4.3U|) . (|4.31f) for F 1 and other well-known objects. 
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6.1 F 1 , isomonodromic tau-function and G-function of Probenius manifolds 

We recall that the genus 1 correction to free energy in topological field theories is given by so-called 
G-function of associated Frobenius manifolds. The G-function is a solution of Getzler equation |3L)| : 
for Frobenius manifolds related to quantum cohomologies, the G-function was intensively studied as 
generating function of elliptic Gromov-Witten invariants (see |24l 131 j for references) . In [23j it was 
found the following formula for G-function of an arbitrary m-dimensional Frobenius manifold: 

G = ln (6.1) 

n m r7 1/48 

llfc=l 'Ikk 

where tj is the Jimbo-Miwa tau-function of Riemann-Hilbert problem associated to a given Frobenius 
manifold ^H]; r]kk are elements of Egoroff-Darboux (pseudo) metric (written in canonical coordinates) 
corresponding to the Frobenius manifold. 

One of the well-studied classes of Frobenius manifolds arises from Hurwitz spaces ^Hj. For these 
Frobenius manifolds the isomonodromic tau-function 77 18 is related to Bergmann tau-function Ty 
(EQgj) as follows |251: 

ri = rj 1/2 , (6.2) 

where / stands for meromorphic function on Riemann surface L. Therefore, the tau-function terms, 
which are the main ingredients of the formulas (|4.3Uj) for F 1 and (|6.1[) for the G-function coincide 
(up to a sign, which is related to the choice of the sign in the exponent in the definition (jl.lj) of 
the free energy). The solution of Fuchsian system corresponding to tau-function 77 is not known 
explicitly. However, the same function 77, being multiplied with certain theta- functional factor, gives 
tau-function of an arbitrary Riemann-Hilbert problem with quasi-permutation monodromy matrices 
which was solved in |26j . 

The metric coefficients of Darboux-Egoroff metric corresponding to Hurwitz Frobenius manifold 
are defined in terms of an "admissible" 1-form (p, defining the Frobenius manifold: 

Vkk = res l^^ • ( 6 - 3 ) 

If, trying to develop an analogy with our formula i|4.3flj) for F , we formally choose <p(P) = dg(P), we 
get rikk = g' 2 (Pk)/2 and the formula (|6.1|) coincides with (|4.3(J|) up to small details like sign, additive 
constant, and the highest coefficient of polynomial V2 arising from requirement of symmetry / <-»• g. 

Therefore, we got complete formal analogy between our expression (|4.3U|) for F 1 and Dubrovin- 
Zhang formula (|6.1|) for G-function. Unfortunately, for the moment this analogy remains only formal, 
since, from the point of view of Dubrovin's theory JH]) the differential dg is not admissible; therefore, 
the metric rj^k = 9 l2 (Pk)/2 built from this differential is not flat, and, strictly speaking, it does 
not correspond to any Frobenius manifold. Therefore, the true origin of the analogy between the 
G-function of Frobenius manifolds and F 1 still has to be explored. 



6.2 F 1 and determinant of Laplace operator 

Existence of close relationship between F 1 and determinant of certain Laplace operator was suggested 
by several authors (see e.g. |27j for hermitial one-matrix model, ^S] for hermitian two-matrix model 
and, finally, |2B] for normal two-matrix model with simply-connected support of eigenvalues, where F 1 
is claimed to coincide with determinant of Laplace operator in the domain with Dirichlet boundary 
conditions). 
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However, in the context of hermitial two-matrix model (as well as in the case of hermitian one 
matrix model |27j ) this relationship is more subtle. 

First, if we don't impose any reality conditions on coefficients of polynomials V\ and V%, function 
F 1 is holomorphic function of our moduli parameters (i.e. coefficients of V±, V2 and filling fractions), 
while detA is always a real function. The Laplace operator A^ which should be playing a role here 
corresponds to the singular metric of infinite volume \df\ 2 . 

This problem disappears if we start from more physical situation, when all these moduli parameters 
are real, as well as the branch points of the Riemann surface L with respect to both projections. In 
this case -F 1 is real itself, as well as determinant of Laplace operator. However, little is known about 
rigorous definition for determinants of such Laplace operators, although such objects were actively used 
by string theorists without rigorous mathematical justification |32l 1341 I31?] . According to empirical 
results of [31], the regularised determinant of Laplace operator A^ is given by the formula 

d6tA/ C\r f f , (6.4) 



AietQB 



where A is a regularised area of L, A^ is Laplace operator defined in singular metric \df(P)\ 2 , B is 
the matrix of 6-periods of L, C is a constant. 

In the "physical" case of real moduli parameters the empirical expression (|6.4|) for lnjdetA^} 
coincides with F 1 (|4,3U|) up to a simple power and additional multipliers. 

Therefore, the relationship between hermitial and normal two-matrix models [2H] on the level of 
F 1 is not as straightforward as on the level of functions F° (F° for hermitian two-matrix mode can 
be obtained from F 1 for normal two- matrix model by a simple analytical continuation |16| I17[ 135] ) . 

From the point of view of determinants of Laplace operators the theorem |2] which tells how the 
Bergmann tau-function depends on the projection choice is nothing but a version of Alvarez-Polyakov 
formula [30]; which describes variation of detA if the metric changes within given conformal class. 



7 From two-matrix to one-matrix model: hyperelliptic curves 

Suppose that di = 1, i.e. polynomial V2 is quadratic. Then integration with respect to Mi in 
can be carried out explicitly, and we get 

Zn = e -^F = c J dMe -NtrV{M) {71) 

where M := Ml, V := V\ and C is a constant. Hence, in this case gives rise to the partition 
function of one-matrix model. 

For di = 1 the function f(P) has two poles of order 1 at 00 ^ and oo 9 ; thus, the spectral curve L is 
hyperelliptic and function f(P) defines two-sheeted covering of CP1. The number of branch points in 
this case equals m\ = 2g + 2; as before, we call them li, . . . , hg+2- The Bergmann tau-function (|3.23|) 
for hyperelliptic curves was computed in |23| : in this case it admits the following, alternative to (|4.6|) . 
(|4.7|) expression: 

T f = A^detA (7.2) 

where 

2g+2 

A:= (h-h), (7.3) 

j<k , j,k=l 
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A is the matrix of a-periods of non-normalized holomorphic differentials on L: 



where 

2g=2 

,2 _ 



k=l 

is the equation of spectral curve L. 

Substituting formula (|7.2|) into ()4.3Uj) . and ignoring coefficient Vd 2 +i (it becomes part of constant 
C), we get the expression 



F^i-lnjA^detA) 12 U ^h)\ 



(7.5) 



which agrees with previously known results E3 El ED] • 
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